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Abstract 

The aspects of phase transitions in the two-dimensional Ising models modified 
by quenched and annealed site disorder are discussed in the framework of fermionic 
approach based on the reformulation of the problem in terms of integrals with 
anticommuting Grassmann variables. 



1 Introduction 

The two-dimensional (2D) Ising model (2DIM) plays important role in the theory of 
phase transitions and critical phenomena due to the analytic results available (in pure 
case) over the whole temperature range [U El 131 Hj- In this report, we review the new 
mathematical methods of analysis and the results so far obtained for 2DIM modified by 
site disorder with application of the anticommuting (Grassmann) integrals. The Ising 
model by itself, in its original formulation, is a lattice model of a ferromagnet presented 
by a set of Ising spins o mn = ±1 interacting with their nearest neighbours along the lattice 
bonds PQ El El II]- The modern approaches to Ising models are merely based however on 
the fermionic path integral reformulation of the problem in terms of the integrals with 
anticommuting Grassmann variables [51 El LTl IHl [9J [10]. The advantage of the use of 
Grassmann variables is that they are canonical variables, as distinct from Ising spins, so 
that one can pass to the momentum space for fermions (5j El [91 [10] . In the pure case, the 
fermionization of 2DIM results the Gaussian fermionic integral for the partition function, 
Z, which in essence means the exact solution of the problem [3, El El [9] . The formulations 
of this kind also admit the interpretation of the 2D Ising model as a lattice quantum field 
theoretical (QFT) problem [Til [121 H31 EH] ■ In particular, the pure 2DIM on a rectangular 
lattice may be presented by the Majorana action with two-component massive fermions on 
a lattice [14J. By doubling the number of fermions, one can pass as well to the Dirac action 
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[H]. The effects of disorder in 2DIM have been extensively studied during last decades 
both theoretically and in the precise Monte-Carlo simulations [H]-[2S]- The disordered 
versions of 2DIM may assume either random modification of the interaction along the 
lattice bonds [TH [121 HI El HI [301 EE], or an admixture of the random nonmagnetic 
impurities at lattice sites pU [191 (221 [251 E2H W\ [28] , also see [321 [331 (Ml ESI [36] . In each 
of this cases one may be interested, motivated by physical considerations, and possible 
applications, in either quenched or annealed versions of disorder. In the quenched version, 
the impurities are assumed to be frozen over the sample. In this case their distribution 
does not depend on temperature and other tuning parameters, like magnetic field, and 
one have to average rather the free energy — (3F = In Z than the partition function Z 
itself, over the impurities [HI [121 DEI E3]- m the annealed version, the impurities may 
be created and annihilated by a variation of external parameters and their concentration 
is governed by the temperature rate and associated chemical potentials [321 E31 [361 EH] • 
In this case, one has to average in Z itself over all states [331 EH ESI [36]. In essence, for 
annealed site disorder, the dilute site can be viewed as being presented by additional (zero) 
component of Ising spin. The resulting model is also known as the spin-1 Ising model, 
or the Blume-Capel model [321 [331 l3B" l 138] . The basic variable in the Blume-Capel model 
is S mn = 0, ±1. The modifications introduced by disorder of any kind typically result in 
the appearance of the additional non-Gaussian terms in fermionic action. Despite of the 
non-Gaussian action in the fermionic integral for Z, the precise results can still be derived 
for disordered 2D Ising models [HI [121 1131 EEH ES] . In what follows, we only consider the 
generic case of the random-site disorder (site dilution) introduced by adding some amount 
of nonmagnetic impurities into a sample, which may be either quenched or annealed, and 
discuss the consequences that can be derived from the fermionic integral representations 
for the partition functions of that models. 

2 The quenched site dilute Ising model 

The basic variable in the pure 2DIM is the dichotomic Ising spin a mn = ±1. The spins 
are disposed at the sites of a regular two-dimensional lattice and interact with nearest 
neighbours along the lattice bonds. The disordered version (quenched site dilution) as- 
sumes that some sites may be nonmagnetic at random. It is suitable to introduce such 
sites by adding the variable y mn = 0, 1 at each mn site, corresponding to the magnetic 
moment of Ising spin at a given site [H]. The resulting hamiltonian is: 

yrnnym+ln@mn@m+ln ~\~ ^2 ymnymn+l@rnn@mn+l\ , (1) 

mn 

where are the ferromagnetic exchange energies; the lattice sites are marked by discrete 
coordinates mn, where m,n = 1, 2, . . . , L, are running in horizontal and vertical directions, 
respectively; we put L 2 — > oo at final stages. For fixed disorder, the partition function 
and free energy are defined by the canonical equations: Z{y} = £ exp(— (3H {y \ a}) = 
exp(—j3F{y}), where the sum is taken over the all possible spin configurations provided 
by u mn = ±1 at each site. The hamiltonian modulus in the Gibbs exponential is: 
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-fiH{y\a} = ^[b 1 ymnym+ln&mn&m+ln ^2 ymnymn+l&mn&mn+l } , (2) 
mn 

where 6^2 = 0Ji,2, and (3 = 1/kT is the inverse temperature. For a typical bond 
weight from Z we write: exp (byy'aa') = cosh(byy') + yy'ua' sinh(6), since aa' = ±1 
and yy' = 0,1. The partition function can then be written in the form: Z {y} = 
R{y} Q {y}, where R{y} is a nonsingular spin-independent prefactor, formed by a prod- 
uct of cosh(byy'), while Q{y} is the reduced partition function: 



Q{y}= s P {II(i )}, (3) 

where ti^ = tanh &x,2, and we assume a properly normalized spin averaging, such that 
Sp (1) = 1 and Sp (a mn ) = at each site. In given case, since we are interesting in 
quenched disorder, we have to average over y mn = 0, 1 rather the free energy —f3\nZ{y} 
than the partition function Z{y} itself. The prefactor R{y} provides only additive non- 
singular contribution like In R{y} to In Z{y} and will be ignored in what follows. The 
problem thus reduces to the averaging of —(3Fq = In Q{y} over y mn = 0, 1 at each 
site. The known device to avoid the averaging of the logarithm is the replica trick: 
[-/3F Q { y }] — [In Q{ y }} = (Q N {y} - l)]jv->o, where [. . .] stands for the average over 
the impurities. In this scheme, one takes iV identical copies of the original partition func- 
tion and average Q N {y} , with formal limit iV — > to be performed at final stages. The 
simplest distribution for the averaging over the impurities is assumed in what follows: 
w {ymn) = pS (1 — y mn ) + (1 — p) 5 (ymn), where p is the probability that any given site, 
chosen at random, is occupied by the normal Ising spin, while 1 — p is the probability 
that the given site is dilute. The averaging of any function like A(y mn ) then results: 
[A(y mn )]=pA(l) + (l-p)A(0). 

The partition function with fixed disorder ([3]) can be transformed into a Gaussian 
fermionic integral following the method of the mirror-ordered factorization for the density 
matrix [H HI [TJ]. Introducing a pair of fermionic (Grassmann) variables a mn ,a mn , we 
write for the horizontal weight: 

1 + tWmnym+ln Cr rnn 'm+ln 

da> mn da mn e mn mn (1 + « mn Z/ mn O" mn ) (1 + t\ a mn ym+ln Cr m+ln) ■ (^) 



1 The situation is different for the annealed case (the Blume-Capel model), where the zero state 
provided by y mn = 0, 1 is rather to be considered as a zero component of the BC spin S mn = 0, ±1 so 
that one have to average over the all three states of S mn = 0, ±1 directly in Z, inside of the logarithm. 
Respectively, the cosine factors (product R{y}) are now to be preserved under the averaging in Z. These 
factors in fact add new degrees of freedom for clusterization of the magnetic sites at low temperatures, 
and are eventually responsible for the appearance of the tricritical point in the Blume-Capel model at 
strong dilution [36 . In fermionic language, the tricritical point is associated with vanishing of the kinetic 
(stiffness) coefficient in the Blume-Capel fermionic action at strong dilution 36J. 
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In a conventional notation, the horizontal weight is now presented as a product of two 
factors, A mn Aj^ +lri , with decoupled spins, taken under the Gaussian averaging. H In a 
similar way, one prepares the factorized vertical weights in the form B mn B* nn+l . At next 
stage, we have to arrange the factors in their global products in order the elimination 
of spin variables be possible [H [91 [H]. The final point is that we have to average over 
a mn = ±1 the product of four factors with the same spin like A^ n Bm n A mn B mn at each 
site, thus passing to a purely fermionic expression. This results the integral [H 

** i-vt n IT] in 



nin 



+ v 2 

' amn 



a 



in^mn + ^2 a m -ln^mn-l + (^l a m-ln + ^2& mn -l) ( a m n + ^mn) }' (5) 



where a mn , a^ n , 6 mn , are Grassmann variables. The disorder parameters y\ in = 0,1 
are still free parameters in the above integral, while Ising degrees being already elimi- 
nated. In turn, integrating out a part of fermionic variables from namely the vari- 
ables a mn ,b mn , we obtain the reduced integral for Q in terms of a* mn ,b* mn [H]. Changing 
notation, a^ n , b* mn — > c mn , —c mn , the integral becomes: 

Q {v} = I IT dc mn dc mn y mn exp ^ ^ y mn c mn c mn + 

mn mn 
+ ( C mn + Cmn) (^l C m-ln — ^^mn-l) ~ Vmn ^2 C m _ ln C mn _ 1 , (6) 

where c mn ,c mn are again Grassmann variables, and we assume: y^ n exp ( y~^c m „c mn ) = 
Umn + CmnCmn, with f/^ n = 0, 1. The integrals ([5]) and ([6]) are still the exact expressions for 
Q{y} originally defined in ([3]). Taking y mn = 1 at all sites, we obtain the 2DIM integral 
for the pure case: 



c c 4- 

I ^mn^mn ' 



Q{ 1 } = / ll dc mndc mn 

mn mn 

( C mn C mn ) (tl c m—ln ^^mn-l) ^1^2 C m—ln^mn—l • (J) 

In particular, the evaluation of the integral (J7J) by transformation to the momentum space 
results the Onsager's expression for Z and InZ of the standard rectangular lattice |14j . 
The advantage of the reduced representation with two variables per site like and (J7J) is 
also that it explicitly illuminates the Majorana-Dirac structures of 2DIM already at the 
lattice level. This can be most easily seen in the pure case. [f| 



2 Let us remember that Grassmann variables (nonquantum fermionic fields) are the purely anticom- 
muting fermionic symbols. Given a set of Grassmann variables a\, 0,2, ... , ajv, we have diCtj +a,jai = 0, and 
a'j = 0. The rules of integration over Grassmann variables (fermionic path integral) were originally intro- 
duced by F.A. Berezin in QFT context [10]. The elementary rules of integration for one variable are [lOj : 
/ dcij ■ a j — 1 , J daj -1=0. In a multidimensional integral, the differential symbols dai, da,2, ■ ■ ■ , dapj are 
again anticommuting with each other and with the variables. Gaussian fermionic integrals are in general 
related to the determinants and Pfaffians. In the field-theoretical language, the fermionic form in the 
exponential under the integral is typically called action. For more comments about Gaussian fermionic 
integrals in a related context also see El [TUJ HH [36] . 

3 In the pure case, the lattice Majorana like action for 2DIM readily follows from (0 by substitution 
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To prepare the iV-replicated integral (EJ), with same set of y mn in each copy, it is 
suitable to multiply first flS]) by factor S(y mn — 1) + S(y mn — 0) = 1, which realizes the 
decomposition over the states y mn = 0,1. We assume that y mn = 1 is realized with 
probability p, and y mn = is realized with probability 1 — p. The averaging of the N- 
replicated integral (jSj) over the disorder within the iV-replica scheme finally results the 
theory with interaction presented by the integral Hi 



TV TV N 

Q N {y}] av = /n n n [p n e5i " + a - *o n 

mn Q =l mn a =l a=l 



?(«) 
mn 



P 



V JV 

nn^^expEiE^ 



(a) 



+ 



mn ct=l 



mn a=l 



1 — p 
P 



A? 



.(o)=(a) -a 



(a) 



u mn ^mn c 



a=l 



where [ s the replicated Gaussian action from (jSJ) for the pure case, see (j7|). The 
effect of dilution is introduced here through the second non-Gaussian term, with bar 
coupling constant go oc The continuum-limit field theory for weak site dilution (RS 
2DIM) that follows from the exact lattice integral flSJ) is commented in more detail in 
[T3] . This corresponds to the low-momenta sector of the exact lattice theory associated 
with (jfJJ) and (jHJ). To extract the effective low-momenta (long- wave) effective action, one 
has to distinguish explicitly the higher and low-momentum lattice fermionic modes in the 
exact lattice action (jSJ). Integrating out the higher-momentum modes in the first order of 
perturbation theory then results the iV-colored Gross-Neveu model (N — > 0) with action 

nr 



S'g-n 
1 



r N 1 

/ d 2 x{ £ [ m N + 5 4 a) (ft + i ft) # 

^ a=l ^ 



v4 o) (-ft + < ft) ^* 

1 — ^1 — ^2 — ^1^2 



(a) 



2 

#iV 



S'iV 



A' 



l-p(B) 1 



2(tit 2 ) c 

y l-p(g) 2 1 
P (A) 2 4(txt 2 ) c 



P < A > /2(tS 



(9) 



where ipi,ip2 are the anticommuting Majorana components, and (7jv are the effective 
mass and charge, respectively. The parameters (A) and (B) are some lattice fermionic 
averages (definite numbers) explicitly calculated in [Tit] . Since the replica limit iV — )■ is 
assumed at final stages, one can put N = and (A) = 1 in mass and charge already in 
(jUJ). The Gaussian part in Qj is the replicated Majorana action, corresponding to the pure 
case, with the mass term modified by disorder. The condition of zero mass will give the 



Cm-in — ► c mn — d m c mn , c mn — > c mn — d n c mn , where d m ,d n are lattice derivatives (momenta). This results 
the Majorana like action S = fhc mn c m n + • • • with mass term and the kinetic part [14] . Evidently, the 
mass parameter will be m = 1 — 1\ — 1% — t\ti. The condition m = defines the critical point of 2DIM in 
pure case [14]. This condition fa = may be rewritten as well in the form: sinh(26i) sinh(2&2) = 1. The 
disordered phase corresponds to positive mass, while the ordered phase corresponds to negative mass. 
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coexisting curve {T c ,p) for quenched site dilution, which is exact for small concentration 
of vacancies, as p — > 1. The analysis of that curve at strong and moderated dilution, that 
is coded in the exact integrals like (0, OH]), has not yet been performed in detail. This will 
probably claim for the advanced methods of approximation like lattice RG or application 
of variational approaches like Hartree-Fock-Bogoliubov. 

The effective continuum-limit N = Gross-Neveu model similar to fl9]), but with 
another itln and gjy, has been originally derived and analyzed by DD-SSL as an effective 
theory near T c for weak bond dilution [TTJ, [T2l [13j E31 [28] - The DD-SSL predictions for weak 
bond dilution, based on the renormalization group (RG) analysis of their effective N = 
Gross-Neveu model in the low momentum sector, with taking also into account some fine 
symmetry effects related to the Kramers- Wannier duality and the CFT interpretations of 
2DIM in the pure case, are the double-logarithmic singularity in the specific heat and the 
logarithmic corrections to the pure-case power laws in other thermodynamic functions, 
as T — > T c [TTj, [121 [131 EU EH]- The derivation of the effective action for site dilution 
in the N = Gross-Neveu form (Q thus supports the idea of the double-logarithmic 
singularity in specific heat and only logarithmic corrections to the pure-case power laws 
in other functions (as in the DD-SSL scheme) also for random-site 2D Ising ferromagnets, 
for weak quenched dilution [TJ]. For more details and a recent discussion of the effects 
of quenched disorder in RB and RS versions of 2DIM along theoretical and experimental 
(Monte-Carlo) lines also see[l9l[22l[23l[2l[25l[lH In conclusion, we 

note that the fermionic integrals like ([5]), ([6]) and ([H]) are still the exact lattice expressions 
for Z (either its reduced version Q). Respectively, one can try other methods of the 
averaging as well as other tools of analysis of fermionic theories of random-site 2DIM 
directly on a lattice, starting from these exact fermionic integrals. 

3 The Blume-Capel model 

The Blume-Capel (BC) model is a classical spin-1 model originally introduced to study 
phase transitions in specific magnetic materials with a possible admixture of non-magnetic 
states. This is a model with the tricritical point at the critical line in (T c , Ao) plane 
[32], [331 EH EH]- The Blume-Capel model can also be viewed as the annealed site-dilute 
version of the ordinary Ising model [36]. The hamiltonian is: 



where S mn = 0, ±1 is the BC spin-1 variable associated with the mn lattice site (m, n = 
1, 2, 3, . . . , L). As distinct from the quenched disorder case, the zero-spin or vacancy state 
of S mn = 0, ±1 is now rather to be considered as a one of the three possible states of spin 
variable [521 E21 ESI ES] • The hamiltonian modulus in the Gibbs exponential is: 



with = /3</i,2 and A = — f3A . The decomposition like S mn = y m n<Jmn, that was 
used in the quenched case, is still possible, but in order to eliminate the Ising degrees by 




(10) 




(11) 



m=l n=l m=ln=l 
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transformation of Z into a fermionic integral it is now more suitable to make use rather of 
the gauge transformation like S mn — > a mn S mn under the averaging, since the product of 
cosines (factor R{y}) is to be included into Z anyhow, before averaging. The Boltzmann 
factors from the Gibbs exponential associated with fTTT]) can be written as (extended) 
polynomials in variables S mn = 0, ±1. The polynomial interpretation is important for 
fermionization [36] . The partition function becomes (with A , = sinh K i , A'j = cosh K j — 

1): 

Z = S P { II II e ASln [(! + A i S mn S m+ m + KS 2 mn S 2 m+ln ) 

^ m=l n=l 

x (1 + A2 S mn S mn +i + A 2 S mn S mn+1 ) j. (12) 

The factorization of local bond weights can again be performed by analogy with the 
Ising case, but now we have to add the even part of the polynomial into the Gaussian 
exponential in the measure [36]. For the horizontal weights, with Grassmann variables 
a m n, a mn , we write: 

1 4~ \\S mn S m+in + i *S' TT1T1 *S' rn ,_|_i rl = 

/ da mn da mn exp{(l -\- X^S mn S m ^ n j & mn (i mn } 
x (1 + a mn S mn ) (1 + Ai a> mn S m +in) j (13) 

and similarly we can factorize the vertical bond Boltzmann weights. The factorization 
( IT3l) makes it possible to pass to the purely fermionic expression for Z in few steps. The 
final fermionic integral for Z appears in the form [36J: 

Z = (2e A cosh K x cosh K 2 ) L / Yl II da mn da mn db mn db mn 

771=1 n=l 

L L 

X exp ^ ^ ^ ] \^Q J mnQ"mn ^mn^mn O'mn^mn 
771=1 n=l 

~\~(t\Clrn-~ln ^2^mn- 1 ) (^tuti ^mn) ^1^2 

9o ^mn^mn^mn^mn GXp ( r )\Q"m—lnQ"rn—ln 72^mn-l^mn-l 

—t\t2 a m ^i n b 

)]}, (14) 
with parameters (where A = — j3A ): 

e- A 



^ = 2cosh^cosh^ 7l = 1 -oo^- = 1 -V 1 -^- ( 15 ) 

The fermionic integral (TH|) is still the exact expression, even for finite lattices, pro- 
vided we assume free boundary conditions both for spins and fermions. The exponential 
in the interaction term can be expanded into a series, which results a finite polyno- 
mial in Grassmann variables. For instance, for particular exponential factor we find: 
exp(— 7 1 a m _ ln a m _i n ) = 1 — 7ia m -i n a rn _i n , and analogously one can expand other factors. 
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The B-C model is thus presented in (fT4l) as a fermionic theory with free-fermion (Gaussian) 
part and a polynomial interaction terms in the action, the highest term in the interaction 
polynomial is of order 8 in fermions. The overall coupling constant in interaction term is 
go oc exp{/3A }, the increasing g Q oc exp{/3A } means increasing dilution. There are also 
additional parameters 7^2 in the interaction polynomial. These parameters come from 
accounting properly the weights related to factors cosh(6iS , mn S , m+ i„) cosh(6 2 5' mri ,S' mn+ i) 
in Z{S} for Blume-Capel, there are no analogs of 71,2 terms in the quenched case. In 
fact, these terms with 71^ are responsible for the existence of the tricritical point in the 
Blume-Capel model. The elimination of the variables like a mn , b mn is now not possible in 
action (Tl^|) . at least straightforwardly, as distinct from the quenched case (cf. (J6])-([7J)), 
just because of the presence of the combinations like a mn a m _i and b mn b mn _i in the 7^2 
terms in the non-Gaussian part of the lattice action (TH1) . 

Despite of the non-Gaussian representation for Z, it is still possible to extract physical 
information by taking the continuous limit (low momenta sector) of the BC lattice action 
like (fbi|) and analyzing it using tools from quantum field theory. The details of construct- 
ing the effective two-component fermionic action at low momenta can be seen in [36]. The 
resulting action includes the Gaussian part, with mass term modified by disorder, and the 
four-fermion interaction of the form (ijjijj | d x i/jd y ip). The condition of the zero effective 
mass already gives the equation for the BC line of phase transitions (critical line) in the 
(T c , Ao) plane, while the effect of the interaction is merely to modify the kinetic terms in 
the action [36]. This also provide grounds to estimate the position of the tricritical point 
[36]. These effects are shortly commented below. 

The critical line is given by the condition of vanishing the mass term in effective BC 
action: rriBc = 1 — ii — £2 — tih +9o = 0, where go is given in (fl5l) . Following [36J, we now 
consider the isotropic lattice case, with tx — t^ — t and K\ = K% = K. The critical line 
is given by msc — 1 + 9o — 2i — t 2 = 0, with t = tanh K and K = (3 J, where = 1/T. 
In notation with K = ^ — ^ and ^-^A , the criticality condition becomes: 

t anh>(i) + 2tanh(i)-l = 5 -i^, (16) 



T 



which may be written as well in the form: 



sinh ( |) = 1 + ~ exp ( ^f) , (17) 



vT/ 2 V T 

which in turn admits the explicit solution for Aq as function of T = T c in the form: 



An = T In 



2sinh(-] -2 



The inverse dependence for T c as function of A can be evaluated numerically by solving 
any of the above equations, which are all equivalent to the condition of the zero mass in the 
effective continuum-limit theory that follow from ( fl^l) . This results the critical line for the 
BC model shown in Fig. 1 in [36]. The critical line is started with maximal T c at the left 
end at A = — oo, which corresponds to the pure case (g = 0), and goes lower as dilution 
increases, with increasing A and go oc exp A °~ 2 . The critical line finally terminates at 
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Chemical potential 




Temperature T C (A ) 




A 


Ref. [33] 


Ref. [35] 


Eqs. (fT6l)-(PT8l) 


0.0 


1.695 


1.714(2) 


1.6740 


0.5 


1.567 


1.584(1) 


1.5427 


1.0 


1.398 


1.413(1) 


1.3695 


1.5 


1.150 


1.155(1) 


1.1162 


1.87 


0.800 


0.800(3) 


0.7712 


1.95 


0.650 


0.651(2) 


0.6135 


1.962 


0.620 


0.619(1) 


0.5776 


1.969 


0.600 


0.596(5) 


0.5531 



Table 1: Numerical values of the critical points (T c (Ao), Ao) in the Blume-Capel model: 
comparison of the results of Monte-Carlo simulations and the equations f ll6p - fll8l) . Note 
that small variation of Ao causes more significant changes in T c (Aq) in the region near 
A = 2, as it is to be expected from ffl6l) - ffT8|) . 

A = 2 at zero temperature. There is no ordered phase at stronger dilution, as it also 
can be deduced from f 1 16 1) - (TT5|) . The theoretical critical line is compared with the results 
of the recent Monte-Carlo simulations for B-C model, see Fig. 1 in [36]. The agreement is 
found to be very good (typically within 1% accuracy) over the whole temperature range 
[36] . The available numerical data for (T c , A ) are also presented (in part) in Table 1. 

The position of the tricritical point at the (T c , Ao) line can as well be estimated from 
the condition of vanishing the kinetic (stiffness) coefficient in the effective B-C action 
associated with ( Tl4"|) [36] . The Hartree-Fock-Bogoliubov method has been applied to 
decouple the four-fermion interaction term in the effective action to extract the corrections 
to the kinetic part [36] • The singular point where the kinetic coefficient vanishes was found 
at (X^AqJ ~ (0.42158,1.9926), in a reasonably good agreement with the results of 
Monte-Carlo simulations for the position of tricritical point: (T t , A 0i t) — (0.610, 1.9655) 
[23], and (T t ,A , t ) ~ (0.609(3), 1.966(2)) [35]. It is in general important that the B- 
C fermionic integral with a non-Gaussian action ffT4"]) finally predicts the existence of 
a tricritical point at the B-C critical line at strong dilution, somewhere close to the 
termination point of that line at Ao = 2. The method of constructing the critical line from 
the condition of zero mass in fermionic integral for Z has been recently extended by Fortin 
and Clusel [37] to the general set of the spin-S' Ising models (S = 1/2, 1, 3/2, 2, 5/2, . . . ) 
[37] . The standard 2DIM and the Blume-Capel models are the first two representatives 
in this set. The agreement of the theoretical predictions for the critical line with the 
available Monte-Carlo data for the spin-S* models was again found to be very good, even 
despite of highly complicated polynomial structures, with many fermions, arising for the 
higher spin- 5* Ising models in the kinetic part of the action [37]. These features may be 
probably understood as an evidence for the well expressed clusterization processes at low 
temperatures in such models, including the generic case of the spin-1 Blume-Capel model. 
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4 Conclusions 



The integrals with anticommuting (Grassmann) variables provide effective tools to analyze 
pure and disordered Ising like spin models in two dimensions. The Ising spin glasses, 
geometry disordered lattices, regularly diluted models, also can be analyzed along these 
lines. In a more general context, it may be noted that there are as well few other important 
physical problems with spins and fermions in two dimensions. The most prominent are 
the quantum Hall effect and the high-T c superconductivity in oxide cuprates and related 
substances. It is interesting that the fine tuning in spin-fermion correspondence as well 
as the effects of disorder seemingly play the role in both cases. It is therefore important 
to understand better the ordering phenomena in terms of fermions in such systems. 
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